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We discuss the recent NA50 J/ if) production data in Pb-Pb collisions, in particular the 
second drop at high transverse energies which correspond to the most central collisions. 
Using a model which relates the J/ip suppression to the local energy density, we show 
that the data can be explained by taking into account transverse energy fluctuations at a 
given impact parameter. Predictions of this model for RHIC are briefly discussed. 
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O ! !• Introduction 

t-H ! 

The rate of J/ip production in p-p, p-A and A—B collisions involving oxygen and 
sulphur projectiles is well understood in terms of a hard production of the cc pair followed 
by nuclear absorption [0,0. In Pb-Pb collisions, however, evidence for an additional 
suppression mechanism, the so-called anomalous suppression, has been obtained by the 
NA50 collaboration |JJ. Furthermore, in the recent NA50 data M a second drop in the 
pattern of the J /if) production occurs at high transverse energy (Et), that is, for the most 
central collisions. The origin of this second drop is the focus of this contribution. 
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2. Improved geometrical model 

It is possible to explain the anomalous J/if> suppression by using a scenario f| which 
relates it to the local energy density e. More precisely, we formulate a simple geometrical 
model in which final state interactions suppress all the J/if>'s originating from cc pairs 
produced in a region where the local energy density e exceeds some critical threshold 
e c . We then assume that the local energy density is proportional to the density n p of 
participant nucleons per unit area (transverse to the collision axis) calculated in a Glauber 
model : the total transverse energy is then proportional to the number of participants, 
as observed experimentally. In such a model, the suppression criterion depends only on 
the impact parameter (through n p ), so that the J/tf> suppression saturates at high E T 
where the impact parameter is essentially zero (see the short-dashed curve in Fig.[l]b). A 
clear deviation from this simple behavior is seen in the most recent NA50 data. 

A better description is obtained by taking into account the fluctuations of Et for a 
given impact parameter b 0. These fluctuations are an essential component of the tail 
of the Et distribution (see Fig.|l|a). As a simple ansatz, we take the distribution of Et at 
fixed b to be a gaussian with mean value (E T )(b) = qN p (h), where N p (b) = J d 2 s n p (s, b) 
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is the total number of participants at impact parameter b and n p (s, b) the corresponding 
density per unit area at the transverse coordinate s. The dispersion of the gaussian is given 
by <j\ = aq 2 N p (b), with a a dimensionless parameter. The values of the fit parameters 
q = 0.274 GeV and a = 1.27 are those determined by the NA50 collaboration [7]. 

Part of the Et fluctuations is phys- 
ical, part is due to the intrinsic res- 
olution of the NA50 electromagnetic 
calorimeter. The latter is given by 

^intrinsic (E T ) / E T = a/y/E^ + P, with 

a ~ 0.2 GeV 1/2 and [3 ~ 0.005 @. 
The physical E T fluctuations <7 p hy S = 

thus differ only 

by 10% from the observed fluctua- 
tions <Je t for central collisions {Et ~ 
100 GeV), so that we neglect the in- 
trinsic resolution in what follows. 

In order to take into account the Et 
fluctuations, we replace the previous 
estimate e oc n(s, b) 0, with the more 
accurate e oc (E T / (E T )(h))n p (s, b) 0. 
Thus for a given impact parameter, the 
energy density is proportional to Et 
[PUfl . With this prescription, the aver- 
age energy density in the collision area 
S, defined by (e) = (1/S) J s e(s,b)d 2 s, 
is proportional to Et/S, in agreement 
with the traditional Bjorken estimate 
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Figure 1: a) Minimum bias ^-distribution, b) 
J/ip over Drell-Yan production ratio as a func- 
tion of Et- Open (resp. closed) circles: NA50 
1996 (resp. 1998) data. Dotted-dashed curve: 
nuclear absorption alone. Full curve (resp. short 
dashes): our model for anomalous suppression 
with (resp. without) Et fluctuations. Long 
dashes: prediction for Au-Au collisions at RHIC 
energies (see text). The vertical line sets the po- 
sition of the knee of the _Er-distribution, where 
the second drop occurs in the ratio. 



3. Comparison with NA50 data 

It is easy to understand that this im- 
provement leads to an increased sup- 
pression at high Et- Indeed, at high 
Et where the geometry of the colli- 
sion is essentially fixed at zero impact 
parameter, the local energy density e 
scales like Et and the region where it 
exceeds e c becomes bigger and bigger 
as Et increases. 

This model reproduces the main fea- 
tures of the J/ij) production pattern observed by NA50; however, it does not provide a 
perfect fit to the data |J. A better description is obtained if we allow for two thresh- 
olds, in line with the idea of successive meltings of the \ an d the J/ip (full curve in 
Fig.|l|b). However we should emphasize that an identically good fit is obtained for a grad- 
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ual suppression above a single threshold. Therefore the structure in the pattern of J /if} 
suppression cannot be interpreted as a signal of the successive meltings of charmonium 
resonances. In any case, we stress that the convolution of the suppression factor with 
the impact parameter distribution (for a given Et) and the energy density profile in the 
interaction region tends to smoothen any threshold, even that associated with the onset 
of the anomalous suppression at lower Et- 

Whatever the chosen scenario and/or the value of the thresholds, we always obtain a 
second drop starting around the knee of the .^-distribution, which directly reflects the 
effect of Et fluctuations. Note that in order to account for the data, all the J/^'s must 
be suppressed at the highest energy densities. 

We have estimated the average transverse momentum squared of the produced J/^'s as 
a function of Et- Analogous computations have been done in Jl2| and we reproduce the 
behavior predicted there: at low Et, {p T ) increases rapidly due to initial state scattering 
of the cc pair |R|; this increase saturates when anomalous suppression sets in, and (p T ) 
eventually decreases above the knee of the Et distribution. Compared to the calculation 
in ]n|, this decrease is amplified by Et fluctuations. This behavior is not compatible 
with the centrality dependence of (p T ) recently measured by NA50 Indeed the data 
show a monotonous increase with Et- Note however that our suppression criterion is 
^-independent, while various arguments lead us to expect that a J/if> with a high pt is 
to be less suppressed. The implementation of such effects is under way [IHJ. 



4. Predictions for RHIC 



The recent PHOBOS measurement |T6[ shows that the total multiplicity is larger by 
70% at RHIC than at SPS in central collisions. We thus assume that the energy density is 
also increased by 70% at RHIC. The above model then yields a parameter free prediction 
for Au-Au collisions at RHIC, which is plotted as the long-dashed curve in Fig.|l]b. In 
order to compare RHIC and SPS results, we have rescaled the RHIC tranverse energy so 
that the knees of the two distributions coincide. In the Glauber model formulas, we use 
the value of the nucleon-nucleon inelastic cross section cjnn = 41 mb at RHIC energies, 
instead of a^N = 32 mb at SPS energies. 

We have assumed so far that the energy density scales with the density of participants, 
i.e. e oc n p (h, s)Et/ (Er)(b); this was a consequence of the observation that multiplicities 
and transverse energies approximately scale with the number of participants at SPS. But 
at RHIC energies, a so-called "hard" component appears which gives rise to an additional 



term, proportional to number of binary collisions |[l^JT7|] . The density per unit transverse 
area of binary collisions is n bin (b, s) oc T A (s)T B (b — s), where Ia(s) is the nucleus profile 
function, and this term must be taken into account in evaluating the energy density. We 
have checked numerically that this does not change qualitatively our predictions for J /if) 
suppression. 

Since the energy density is significantly larger at RHIC than at SPS, the anomalous 
suppression sets in earlier. Around the knee, the suppression for Au-Au collisions is about 
2.5 times greater than at SPS for Pb-Pb collisions, so that the effect of Et fluctuations, 
although still visible, is much reduced. In a smaller system, e.g. Ca-Ca, the J /if) suppres- 
sion would be roughly the same at RHIC as the Pb-Pb system at SPS. Note, however, 
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that the present model ignores the possibility that several cc pairs may be produced in 
the collision. As recently shown |T3 , this may lead to an enhancement of J / ip production 



at RHIC, and may mask the suppression mechanism discussed here. 
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